Classical Casimir interaction in the plane-sphere geometry 
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We study the Casimir interaction in the plane-sphere geometry in the classical limit of high tem- 
peratures. In this limit, the finite conductivity of the metallic plates needs to be taken into account. 
For the Drude model, the classical Casimir interaction is nevertheless found to be independent of the 
conductivity so that it can be described by a single universal function depending only on the aspect 
ratio X — L/R where L is the interplate distance and R the sphere radius. This universal function 
differs from the one found for perfect reflectors and is in principle amenable to experimental tests. 
The asymptotic approach of the exact result to the Proximity Force Approximation appears to be 
well fitted by polynomial expansions in In a;. 
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The Casimir effect arises due to the confinement of the 
quantum fluctuations of electromagnetic fields between 
reflecting bodies. Its discussion is often focused on the 
ideal case introduced by H. Casimir |T| of a pair of per- 
fectly reflecting parallel plates at zero temperature. In 
this case, the energy is given by a universal expression 
8q — ~'hcTT^A/720L^ depending only on geometrical pa- 
rameters, the mirrors' separation L and surface A, and 
fundamental constants, H and c. However this idealiza- 
tion does not hold for any experiment and reliable de- 
scriptions of the Casimir interaction have to account for 
the optical response of matter [2 12] . 

In the present letter, we want to study another case 
where the Casimir interaction becomes universal, namely 
the high-temperature limit indicated in the following by 
a subscript HT. This limit has been thoroughly studied 
for two perfectly reflecting parallel plates P'J5| where the 
free energy J^^rp = —({3)kBTA/8TTL'^ depends only on 
temperature T and geometrical parameters L and A {C, is 
the Riemann function and C(3) ~ 1.202; the superscript 
P refers to perfect mirrors). When the Drude model is 
used for describing the effect of conduction electrons in 
metallic mirrors, the free energy J-^^p = T^r^/2 is one half 
of that obtained for perfect mirrors |7j-i9j (the superscript 
D refers to the Drude model). 

This reduction by a factor 2 at large temperatures, 
independent of the parameters of the Drude model (see 
below), is confirmed by microscopic descriptions of the 
interaction between two metallic bulks p!0Hl2] . This dif- 
ference can in principle be tested at room temperature 
by measuring the force at large distances. Experiments 
are difficult in this domain because the force decreases 
with distance. However, experiments have recently been 
performed at distances up to 7/im where the thermal ef- 
fect is large. The results have been interpreted by the 
authors |13j as being in agreement with the Drude model 



prediction once an electrostatic patch contribution is sub- 
tracted [M] . This conclusion stands in contradiction with 
the results of other Casimir force measurements, per- 
formed at smaller distances, up to 0.75/im, which were 
interpreted by their authors as excluding the dissipative 
Drude model and agreeing with the lossless plasma model 
[15] . This contradiction remains a matter of debate (see 
discussions and references in [El Hi]). 

It is also worth stressing that most precise experiments 
are performed in the plane-sphere geometry. The evalua- 
tion of the force is often done through the so-called Prox- 
imity Force Approximation (PFA) [18] which amounts to 
averaging the force between parallel plates over the dis- 
tribution of local inter-plate distances. This trivial treat- 
ment of geometry cannot reproduce the rich relation ex- 
pected between Casimir effect and geometry [TO ] 120 ] [22 ] . 
Theoretical treatments going beyond the PFA have re- 
cently been proposed, in particular through a multipolar 
expansion well adapted to the plane-sphere geometry. In 
the following, we will use results known for perfect or 
metallic mirrors coupled to electromagnetic fields [23] [24]. 

The aim of the present letter is to investigate the clas- 
sical Casimir interaction for the geometry of a plane and 
a sphere with arbitrary aspect ratio. We will focus our 
attention on the models of perfect mirrors and dissipa- 
tive mirrors and use scattering methods discussed in [23] . 
The results will turn out to be given by universal func- 
tions of the aspect ratio, for dissipative as well as perfect 
mirrors, with the two functions differing from each other. 
The limit of high temperatures will allow us to reach a 
much better numerical accuracy than in previous works 
[24) . This improvement in the accuracy will be used to 
sharpen the discussion of the asymptotic approach of the 
exact result to the PFA. This asymptotic behavior will 
appear to be well fitted by polynomial expansions of the 
logarithm of the aspect ratio. 



We start from the Matsubara formula giving the 
Casimir free energy T between a plane and a sphere |24j 



T ^ fcBT^lndetX>(i^„ 



V = X -Tl2er^'^Tlxe-^'^ 



(1) 



Here ^„ = l-Knk^T jh are the Matsubara frequencies, and 
the prime indicates a sum over integers n = . . . oo with 
a factor ^ for n = 0; TZx and 7?.2 are the reflection oper- 
ators on the plane and the sphere respectively, discussed 
below in terms of Fresnel reflection amplitudes and Mie 
scattering amplitudes; e~^^ describes the one-way prop- 
agation on the distance C = L-\- R between the center of 
the sphere and the plane with R the radius of the sphere 
and L the minimal distance between the plane and the 
sphere; K, = -^/fc^ -|- ^2/^2 jg ^j^g longitudinal wavevector 
after a Wick rotation from real to imaginary frequencies, 
with k the modulus of the transverse wavevector. 

The nonzero Matsubara frequencies ^n^o ^ ink^T /h 
increase with temperature and their contributions to the 
sum ([l]) become exponentially small because of the prop- 
agation factor e~'^'-' with JC > ^/c > 2'Kk-QT / ?ic. At the 
high temperature limit k^T 3> he/ C, the free energy is 
dominated by the first term n = in the sum ([I]) and, 
therefore, it is proportional to temperature 
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•D(O) = Z - ■Rs{Q)e'^^np{Q)e-''^ 



$ = — lndctX>(0) 



(2) 



$ is a dimensionless function of the parameters which 
does not depend on temperature. Hence the entropy S is 
independent of temperature and the energy E vanishes 






(3) 



This implies that the classical limit of the Casimir inter- 
action has a purely entropic origin [25l [26j . 

The Drude model is the simplest description of a dis- 
sipative mirror. It corresponds to a dielectric function 
e [i^] = 1-f tr/^ at imaginary frequencies uj = i^, without a 
magnetic response (/i = 1). a [i^] = u!p/{'j + £,) measures 
the conductivity associated with conduction electrons 
(the SI conductivity is eocr). The Drude model is char- 
acterized by two dimensional parameters, the squared 
plasma frequency Wp proportional to the density of the 
conduction electrons, and their relaxation rate 7. It 
meets the important property of ordinary metals to have 
a finite static conductivity (Tq — ^nll- The lossless 
plasma model 7 — >■ corresponds to an infinite value 
for (Tg, in contradiction with tabulated optical data at 
low frequencies [Ul |25| . 

The Fresnel reflection amplitudes tte and ttm for 
plane waves with polarization TE or TM are well known 



from studies of the plane-plane geometry [7H5] . They are 
Ti^g — > — 1 and r^yi — > 1 for perfect mirrors whereas they 
go to rJj^B — >• and r^y^ -^ 1 for the Drude model at 
low frequencies. For the lossless plasma model, they are 
found to interpolate between perfect and Drude models, 
a property found also below for the Mie amplitudes. The 
disappearance of the TE reflection amplitude is responsi- 
ble for the ratio 2 between the classical limits of the free 
energies obtained from perfect and Drude models in the 
plane-plane geometry. 

A similar discrepancy arises also for the Mie ampli- 
tudes ai and bi which describe the scattering on the 
sphere |24j . For perfect mirrors and low frequencies, we 

P (-1)^(£+1)|(2<+1) p a^l 

get Qf, — g(2i+i)\\(2i—i)\i. ana o^ ~ ~£^i- -Lne electric 
and magnetic multipoles ai and bg progressively decrease 
with increasing multipole index £, but they have the same 
power law in the reduced frequency ^ — S,R/c for a given 
£. The situation is different.for the Drude model where 
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((To ~ <7f)R/c is the re- 



duced frequency associated with the static conductivity) . 
Incidentally, the results for the lossless plasma model are 
found to interpolate between perfect and Drude models: 
bg has the same form as b^ (resp. b^) for large spheres 
(resp. small spheres). The crossover between these two 
regimes is determined by the parameter ujp being respec- 
tively larger or smaller than 1 (wp = u-pR/c is the reduced 
plasma frequency). 

The function $^ is a universal expression of the aspect 
ratio for perfect mirrors, as there exist no other param- 
eters. The discussion of the preceding paragraph entails 
that it is also the case for the function ^^ calculated for 
Drude mirrors [21J , although the latter is associated with 
dimensional parameters like ctq. As a matter of fact, bg 
is negligible with respect to ag for any given value of ^, 
and it follows that the classical high temperature limit is 
the same as if b^ would simply be zero. In spite of the 
fact that ^^ and $^ are universal functions of the aspect 
ratio, these two functions differ. As another important 
consequence of this discussion, the lossless plasma model 
result does not obey the universality property met by 
perfect as well as Drude mirrors. For this reason, we will 
not discuss this model any longer in this letter. 

In the remainder of this letter, we present the results 
of the numerical evaluation of the functions ^^--^ which 
determine the classical Casimir free energy between a 
plane and a sphere, for Drude and perfect mirrors respec- 
tively. We write these functions through a comparison 
with the often used PFA approximation (the superscript 
D, P means D or P) 
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The first factors in the expressions of $^'^ are the PFA 



results while the second ones, g^'^, represent deviations 
from PFA. The factors g^'^ go to 1 for large spheres x — >■ 
0, which confirms the validity of PFA as an asymptotic 
approximation. The fact that ^^-^ and g^'^ depend 
on a single geometrical parameter, the aspect ratio x, 
corresponds to the universality discussed above. 

The explicit numerical evaluation proceeds as in |24j . 
Due to the relatively simple form of the matrices in the 
low frequency limit, the numerics can be pushed farther 
than for the general case. Precisely, the computation can 
be pushed to a larger maximum value €niax for the mul- 
tipole index £, and it follows that the numerical results 
are now more accurate for small values of the aspect ra- 
tio X. The results of the numerical evaluations shown on 
Fig. [T] have been calculated with ima.^ up to 5000 at our 
minimal value for the aspect ratio 10"'^. The numerical 
accuracy 6g is estimated to be of the order of lO^** at 
X = 10^3 and of 10"* at a; = 2 x 10"^. 
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FIG. 1: Classical (high temperature) Casimir interaction be- 
tween plane and spherical mirrors, represented by the ratio 
Q^'^ of exact result to PFA result versus the aspect ratio 
X = L/R; the full red line corresponds to Drude mirrors and 
the dashed blue line to perfect mirrors [Colors online]. 

Though both of them are universal, the functions 
g^{x) and g^{x) drawn on Fig. Ill do not coincide. In 
other words, the ratio $^(x)/$°(x) = 2g^{x)/g^{x) be- 
tween Casimir free energies calculated for perfect and 
Drude mirrors is itself a function of x. This ratio reaches 
the value 2 at the PFA limit (a; — > 0) , as known from the 
studies on plane-plane geometry j7j, and the value 3/2 at 
large distances {x — >• oo) fM\ . 

We finally discuss the asymptotic approach to PFA for 
small values of x. To this aim, we introduce another 
representation of the deviation from PFA 

gD.P(a;) = l-Ha;/?D'P(x) . (5) 

/3°'^(a;) is the slope of the line which joins the points 
(0, 1) to (a:, g^'^{x)) on the plots of Fig. [l] It may also 
be thought of as an additive correction to PFA in Q 

$D'P(x) = CD'P('^+/3D.P(x)') . (6) 

The numerical accuracy 6g discussed above is now trans- 
lated into a numerical accuracy 6/3 = Sg/x for /3. This 



accuracy Sf3 is estimated to be of the order of 10 ^ at 
X = 10-3 and of 5 x IQ-^ at a; = 2 x lO'^. The follow- 
ing discussions of the asymptotic approach to PFA are 
possible only thanks to this good accuracy. 

The quantities f3^'^ are shown as the upper plot on 
Fig. |2]for Drude and perfect mirrors, as functions of In a;. 
The lower plot on Fig. 2 represents the quantities — ^-^^^^ 
which appear in the Casimir force 
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FIG. 2: Deviation from PFA of the classical Casimir inter- 
action between plane and spherical mirrors, described by the 
additive correction fP'^ (x) versus In a; (upper plot) and the 



d^S"- 



a: j„ (lower plot); the 



losarithmic derivatives . , — .x. . 

^ a In x ax 

full red lines correspond to Drude mirrors and the dashed blue 
lines to perfect mirrors [Colors online]. 

The plots on Fig. [2] suggest that the asymptotic ap- 
proach to PFA is well described for small values of x by 
a polynomial expansion in the variable In a;. In order to 
check this idea, we have performed best fits of the nu- 
merical results with the following trial functions, defined 
with the same number of parameters, 

/3i (a:) = ao + fli In a; + 02 In .a; + 03 In x , 



jiySx) = 60 + ^\x ^' b2X^ + h'iX^ , 
Pm{x) = Co + Ci Ina; + C2a; + c^x^ 



(8) 



The first one, /3i, is a polynomial form in the variable In a; 
as suggested by Fig. |2] (see also [29]). The second one, /3x, 
is a polynomial form in the variable x, which was suffi- 
cient to represent functions p obtained from older calcu- 
lations with a lesser accuracy [321 130] ■ It corresponds to 



the free energy being a Laurent series in the variable x. 
The third one, /3ni, is a mixed form which results, in the 
classical limit (high temperatures), from the assumption 
[5T] that the force (not the free energy) may be written 
as a Laurent series in x. 
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FIG. 3: Best fits of the numerical results with the trial func- 
tions ([8| for the Drude (red lines, upper plot) and perfect 
mirrors (blue lines, lower plot); the full, dotted and dashed 
lines correspond respectively to /3i, /3x and fim', the full circles 
correspond to numerical points not used for the fits; the insets 
zoom on the smallest values of x [Colors online]. 



In each case, the best fit functions are defined from the 
numerical results obtained in the window Inx G [—5, —3]. 
Their comparison on Fig. |3] with the numerical values 
obtained in a broader window shows that /3i is a better 
representation of the numerical results than /3m, itself 
better than /3x- In particular, the insets on Fig. [3] show 
the deviations of the best fits of /3^ and /3m from the 
exact results at the smallest values of the aspect ratio 
lux £ [—6.2,-5]. The assumptions that the free energy 
or force be a Laurent series in x cannot be considered as 
generally valid. 

The difference between the cases of Drude and per- 
fect mirrors is in principle amenable to experimental 
tests. Direct comparison with exact calculations would 
require measurements at large distances (high tempera- 
ture limit). This challenge could be met with experimen- 
tal parameters not so far from those in the experiment 
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